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Abstract: Using a combination of Kawashima- and Goodman-type energy esti- 
mates, we establish spectral stability of general small-amplitude relaxation shocks 
of symmetric dissipative systems. This extends previous results obtained by Plaza 
and Zumbrun [9] by singular perturbation techniques under an additional techni- 
cal assumption, namely, that the background equation be noncharacteristic with 
respect to the shock. 



I. Introduction 
Let us consider the one-dimensional hyperbolic system with relaxation 

(1) w t + F{w) x = Q{w) 

for the unknown w = w(x,t) G R N , x G R,t > 0. Here F G C 2 (R N ; R N ) is such 
that dF(w) has N real distinct eigenvalues for any state w under consideration and 
Q G C X (R N ; R N ) has the structure Q(w) = (0 n ,q(w)) where q G C l {R N ; R r ), 
r = N — n. Additionally, we assume the function q to have a relaxation structure: 
let w = (u, v) G U x V C R n x W , 

i. there exists a C 1 function v * : U C R n — > !R r such that q(w) = in U x V if and 
only if v = v*(u) where w = (u, v) G M n x R r ; 

ii. for any u &U, all of the eigenvalues of d u q(u, v*(u)) have negative real part. 

As a consequence, it is natural to introduce the corresponding relaxed hyperbolic 
system of conservation laws, formally obtained by considering the first n equations 
of ([1]) and substituting the variable v with the equilibrium v = v*(u) 

(2) u t + f*(u) x = where /»:=/(«,«»). 



System (TjQ) possesses smooth traveling wave solutions corresponding to shock waves 
of the relaxed system (121) at least in the small-amplitude casej* Existence of such 
special solutions has been given for specific models in the large-amplitude case - 
for example, the Broadwell model, [1] — , or for general relaxation system in the 
small-amplitude case, see [TT1 17]. By changing frame, such travelling wave can be 
assumed, without loss of generality, as stationary solution of (JT]), i.e. solution of the 
form 

(3) W = W(x), W(±oo) = W±. 

where W± = (u±, v*(u±)) with u± denoting the state connected by the corresponding 
relaxed shock wave. 

The next natural question to answer is whether such steady states are stable or 
unstable. The equation for the perturbation w := w — W is 

w t + (F(W + w)- F(W)) x = Q(W + w)- Q(W). 

and the corresponding linearized equation is 

(4) w t = Cw:=- (dF(W)w)' + dQ(W)w. 
Thus, the linearized eigenvalue equation is 

(5) (A / - C)w = A w + (dF(W)w)' - dQ(W)w = 0. 

From now on, we consider equation ([5j) in the Sobolev space H 1 and we say that 
A G C is an eigenvalue of C if there exists a function w G H 1 \ {0}, such that ([5]) 
holdsfj 

By differentiating the equation satisfied by the profile W, we get 

(dF(W)W')' - dQ(W)W = 0. 

In the noncharacteristic case, i.e. dF(W±) invertible, W decays exponentially fast 
to zero as \x\ — ► oo. Thus, A = is an eigenvalue of the linearized operator L. Hence, 
instability is related to the presence of non-zero eigenvalues with non-negative real 
part. 



1 In all of the paper, we will not consider the case of nonsmooth relaxation shock profiles, i.e. 
exhibiting subshocks. This is not restrictive since the smallness assumption of the profile, usually, 
guarantees also smoothness. 

2 On the region Re A > and A ^ that we will consider, and under our hypotheses (A1)-(A2) 
below, H 1 spectrum agrees with L p spectrum for any 1 < p < oo; see [5]. 
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Definition 1.1. The stationary solution W is spectrally stable in H 1 , if for any A G 
C \ {0}, Re A > 0, whenever w G H 1 solves the resolvent equation (jSJ) then w = 0. 

Results in [7J[8] show that, under additional assumptions, spectral stability implies 
both linear and nonlinear orbital stability. Hence, determining whether spectral 
stability holds or not is the key issue for determining nonlinear stability of relaxation 
shock profiles. The aim of the present paper is to prove a general result on spectral 
stability of relaxation shock profiles assuming smallness of e := \ W+ — W-\. 

We divide the assumptions into three groups: A, B and C, referring, respectively, 
to the relaxation system (JTJ), to the the relaxed system ([2]), to the relaxation shock 
profile ([3]). These are imposed on a small neighborhood WcR" about an equilib- 
rium point wo = (uo, vq), q(wo) = 0, or equivalently Vq = v*(uq), and a neighborhood 
U* of uo such that the graph of v* over U* is contained in W. 

Al. There exists a smooth function A = A (W) from W C M. n to the set of real 
symmetric positive dehnite matrices such that (A dF)(W) is symmetric for any W 
under consideration, and 

(6) Re(w,(A dQ)(W±)w) < -c\U ± w\\, 
for some c > 0, where U± w := dQ(W±) w. 

A2. (Shizuta-Kawashima condition) There exists a smooth skew-symmetric matrix- 
valued function K, depending on dF, dQ and A such that 

(7) Re (KdF - A dQ) > 0. 

Bl. There exists a smooth function cto = clq(u) from U* C M p to the set of real 
symmetric positive definite matrices such that aodf*(u) is symmetric for any u 
under consideration and a b is positive semidefinite, Re a &o > 0, where b as 
defined in (jHJ) is the associated Chapman-Enskog viscosity. 

B2. (reduced Shizuta-Kawashima condition) There exists a smooth skew-symmetric 
matrix-valued function k, depending on df*,b and a such that 

(8) Re (kdf* - a b ) > 0. 

B3. (Simplicity, genuine nonlinearity of principle equilibrium characteristic) There 
exists Co > such that there is a single eigenvalue ao(u) of df*(u) that has absolute 
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value < Co on W , with all others of absolute value > 2cq. Moreover, 

(9) da (» • r Q (u) =: r](u) ^ 

on 14* , where ro(u) denotes the unit right eigenvector associated with ceo(u). 
C. There exists C > such that for any x G M there hold 

(10) |W"| i00 < C\W + - W_| 2 , \W"(x)\<C\W+-W-\\W'(x)\, 



and 



w 

— — + sgn (r))r 



< C\W+ - 



Remark 1.2. The apparently restrictive A1-A2, B1-B2 in fact all follow from the 
standard assumptions that (i) there exist a positive definite symmetrizer Aq, A$dF 
symmetric, that at equilibrium points simultaneously symmetrizes dQ, A dQ sym- 
metric (weak simultaneous symmetrizability) , and (ii) at equilibrium points, no 
eigenvector of dF is in the kernel of dQ (genuine coupling); see Lemma lA.ll Ap- 
pendix El These two assumptions hold quite generally in applications, in particular 
for discrete kinetic equations and moment closure systems |10j . Assumption B3 is 
standard and easily checked. 

Remark 1.3. Assumption C is satisfied for a family of profiles near wq if: 

(i) dF is invertible and a (u ) = (see Appendix of [7j). 

(ii) dF is constant, a (wo) — 0, and dimension bounded, e.g. in the case of 
discrete kinetic models with upper bound on the number of modes (see [8]). It has 
been shown to hold also for the infinite-dimensional case of the Boltzmann equation 

®- 



With these assumptions, our main result is as follows. 

Theorem 1.4 (Spectral stability). Under assumptions A1-A2, B1-B2-B3, and C, 

for e := \W + — W-\ sufficiently small, the relaxation shock W is spectrally stable. 

As the argument is somewhat complicated, it may be helpful to outline here the 
structure of the proof. We start by carrying out the following by-now-standard 
"Kawashima-type" energy estimates on the relaxation system. 
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Proposition 1.5. Assume hypothesis A1-A2 and C. Let A G C such that Re A > 
and let w be a solution of Then for e sufficiently small, there hold: 

(12) Re A \w\ 2 2 + |nw|^ 2 + \w'\ 2 2 < Ce 2 \u\ 2 2 (Kawashima estimate) 

(13) ReA<Ce 2 , |ImA|<C£. 
where II := II + + II_ . 

Evidently, it remains only to obtain estimates on the equilibrium variable |w| 2 2 . To 
this end, we carry out an approximate Chapman-Enskog expansion, keeping track 
of error terms, to obtain an effective viscous system for u of the same symmetric 
dissipative type, but with error terms coming from higher derivatives. Applying 
Goodman-type energy estimates to the integrated version of this reduced system, 
following [5], we obtain the desired bounds on \u\ 2 modulo errors consisting of 
higher-derivative and dissipative terms (denoted by v in Section [3]). Observing that 
these, by ( fl2|) . may be absorbed in lower-order and equilibrium terms, we are done. 

More precisely, we establish the following bounds on the reduced system. Here 
we use the following notation for the Vy'-weighted L 2 -norm 

1/2 

i /TT7777 i / / i i9 n„ii / \ r \ 

\Z\ 



\VW r \A L z = U\z\ 2 \W'\(x)dx) 



The space of functions with bounded \ ■ \ w , will be denoted by L 2 r Since W is 
bounded, there holds \z\„ r . < C\z\ . for some C > (a natural choice is C : = 
\W'\^). Hence L 2 C L 2 t with continuous injection. The opposite inequality is 
false since W decays (exponentially fast) to zero as \x\ — > oo. 

Proposition 1.6. Assume hypothesis B1-B2 and C. Let z be defined as follows 

z{x) := / u(y)dy. 



Then, for e sufficiently small, there holds 

(14) ReAlzl 2 , + Id 2 , < C\z\ 2 . 

\ I L 2 L 2 W 

Thanks to (114j) , it is sufficient to control the weighted norm \z\ ,. 

Proposition 1.7. Under assumptions B1-B2-B3 and C, for e sufficiently small, 
there holds: 

(15) ReA|z| 2 2 + \z\ 2 < C(e\u\ 2 l2 + e\TIw\ 2 l2 + e~ 1 \w'\ 2 l2 ) (Goodman estimate) 
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where II := II + + II_ 



Theorem 11.41 is an immediate consequence of Propositions 11.51 11.61 and 11.71 



Proof of Theorem Combining ffTJZj) and (fT5|) . we obtain 



With (HU), this gives 



kl 2 <ce|u| 2 ,. 



|«P <C|z| 2 <Ce 2 \u\\ . 



showing that, if £ = \W + — W_\ is small enough, w = 0. □ 

Discussion and open problems. We remark briefly on the setting of these re- 
sults. Small-amplitude existence and stability were shown in [TT], [7] and [9] under 
the additional noncharacteristicity assumption det dF ^ 0, or, equivalently, the con- 
dition that characteristic speeds of the background system do not vanish relative to 
the shock speed. This hypothesis suffices to treat simple model problems such as 
the Broadwell or Jin-Xin equations. However, as discussed in [8], it is unrealistic 
for models derived by discretization or moment closure from kinetic equations, since 
these may possess characteristics of any speed. Thus, it is highly desirable to re- 
move this technical hypothesis, as we do here. The combination of Goodman- and 
Kawashima-type energy estimates was used in [5] to treat stability of viscous shock 
profiles for systems with real viscosity. A similar, but more complicated argument 
combining these ingredients was used in [6] to treat stability of Boltzmann profiles. 
These results motivate the present analysis, which essentially interpolates between 
the two. 

Interesting open problems are verification of linearized and nonlinear stability in 
the same setting, assuming spectral stability, and the direct verification of C using 
stability estimates together with known bounds on the profile for the reduced system 
following the philosophy set out in [6]. 

Notations and (very) basic tools. Given ti>i,ii>2 : K — ► C n , we denote by 
(wi,W2) the scalar product defined as follows 



(wi,w 2 ) := / w\(x) ■ W2(x) dx 



where w denotes the complex conjugate vector of w. Given A, n x n matrix with 
complex entries, there holds 

Re (w, Aw) = - (jw, Aw) + (w, Aw)^j = (w, A*w) 

where A* := (A + A*)/2. 

If S : R — > W axn is such that S(x) is symmetric for any x, then 

Sw') = l w ■ Sw' dx = — / Sw' -wdx— / w ■ S' w dx = —(w, S w') — (w, S' w) 
Jm Jr Jr 

Hence 

(16) Re(w,Sw') = --(w,S'w). 
Similarly, if K is skew-symmetric, then 

(17) Im(w',Kw) = ~(w,K'w) 

In particular, if K is constant, (it/, Xiy) is a real number. 

From here on, we will denote with 0(1) any function of x, W and A, locally 
bounded in {(x, W, A) : Re A > 0}. As a consequence, given the functions f,g G L 2 
and h G L 2 , the following estimates hold 

(is) |(/,o(i)r ? )|<qru (vlf^+v- 1 ^) 

(19) | (n, 0(1) W'g) | < O (n|/f w , + iTWV |< 2 ) 

where 77 is any strictly positive constant and O is a constant independent on 77. 

2. Estimates for the full system 

Lemma 2.1. Lei e := |W+ — W_| and assume hypothesis Al and C. Let A G C snca 
i/iai Re A > let w be a solution of (TJP and Zet be any constant skew- symmetric 
matrix. Then for e, rj > both sufficiently small, there hold 

(20) Re X\w\ 2 + \Uw\ 2 < Ce 2 \w\ 2 o , 

(21) ReA|w'|2 2 - Re (it/, A dQ w') < C e 2 \w\ 2 h1 ; 

(22) Re (w',KdFw') < C (s 2 \w\ 2 l2 + n" 1 \Uw\ 2 l2 + (e 2 + v)W\ 2 l2 ) 
where IT := Il + + Il_ and C denotes a constant independent on e and rj. 
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Proof. Taking the scalar product of A (W)w against (jSJ), we obtain 

(23) X(A w,w) + (A w, (dFw)') - (A w, dQ w) = 0. 
Hence, using (fl6l) . we get 

Re\(A w,w)-Re(w,A dQw) < -Re (A w, d 2 F W'w) + ^Re (w, d(A dF) W' 
Since A is positive definite, there holds for some C > 

(24) ReAMj 2 -Re(w,A dQw) < C \W'\ LOO \w\\ 2 . 



Let us set 



\W - WJ \W - W I 

Then there holds, for some C > 0, 

|$(W) - (A Q dQ)(W)\ < C\W -W-\\W -W+\. 

Therefore 

Re(w,A dQw) < Re(w,$(W)w) + C\W -W-\ LOO \W - W+\ LOO \w 
Thanks to (jSJ), we get 

Re( W ,A dg W ) < -c\Uw\ 2 l2 +C\W-W.\ LOO \W-W + \ Lao \w\ 2 L2 
for some C, c > 0. Hence, using ( 1241) . we obtain 



2 

Z.2 



(25) ReAM 2 L2 + |nH* a < C (\W - W + \ LOO \W - W.\ LOO + \W'\ LOO ) \w\\ 2 . 

In term of e, we get the 0-th order Friedrichs estimate ( 1201) . 
Differentiating (jHJ) with respect to x, 

(26) A w + (dF(W»" - (dQ(H/)w)' = 0. 
Taking the scalar product of Aq(W) w' against fl26l) . we get 

(27) A(A - (w',A dQw') = -(A w', (dFw)") + (A w', d 2 Q W'w). 
Since 

(A w', (dFw)"} = (A Q w\d*FW'W w) + (A w', d 2 F W" w) 
+2(A w',d 2 FW'w'} + (w',A dFw"), 
taking the real part and using (TIB"]) , we obtain 

Re(A w', (dFw)"} = (w',0(l)W'w) + (w',0(l)W"w) 



+(w',0(l)W'w') - ~(u/,d(A dF) W'w'). 
Hence, the following estimates holds 

\Re(A w', {dFw)")\ <C(\W'\ LOB + \W"\ LOa ) 

\Re(A w',d 2 QW'w)\ < C \W'\ LOO \w\* . 
Therefore, from fl27|) . using ffTUj) . we deduce 

ReA|w'|* 2 -Re«A dQw') <C|W*| i00 \w\ 2 b1 . 

Thus in term of s, we obtain the 1st order Friedrichs estimate ( |2T|) . 

Now, let K be any constant skew-symmetric matrix. Applying K to the resolvent 
equation ([5]) and multiplying by w', we get 

\(w', Kw) + (u/,fr(dFw)') - (w',KdQw) = 0. 

Taking the real parts and rearranging the terms, we obtain 

Re (w',KdFw') = -Re (A (w',Kw)) - Re (w', K d 2 F W' w) + Re (w',KdQw). 

Hence, thanks to ( JT71) . there holds Im (it/, Kw) = 0, since K is constant. Therefore, 
for Re A > 0, we obtain 



Let us set 



| Re (A (w', Kw)) | = Re A |Re (w', Kw) \ < C Re A \w\ 2 x . 



\W - W.A \W - W I 



Then there holds, for some C > 0, 

-lfdQ(W)| < C7|W- W+\. 

Therefore 

Re (u/, KdQw) < Re (w', to) + C\W - W_| iOQ | W - W + | i00 
Thanks to (jBJ), we get 
Re(w',tfdQw) < C^in^l^ +77|w'|2J +C|H/-H/_| i0O |H/-H/ + | iOO |w|2 2 

where rj > is a positive constant to be chosen later on (small enough). Hence, the 
following estimate holds 

Re (w',KdFw') = C {(ReX + \W'\ LOC + \W -W-\ LOO \W -W+\ loo )\w\ 2 h1 

+ r]~ l \H.w\ 2 + r]\ w '\ 2 o] 
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By (120|) . ReA < Ce 2 , hence, in term of e, we get ( 1221) . This concludes the proof of 
Lemma 12.11 □ 

Assuming, in addition, hypothesis A2, we prove estimates ffT2l) and (TTB"]) . 

Proof of Proposition \1.5[ Thanks to the small amplitude assumption, it is possible 
to choose K = K(W + ) constant in the Shizuta-Kawashima condition (J7J), since 
this is an open condition so persists under small perturbations. Hence, summing 
estimates dSUD-dH]) with (J22]), we obtain 

Re A \w\\ + \Uw\\ + \w'\\ <Cs 2 \w\ 2 l + C (1 + rT VMl + C(e 2 + r l )\w'\ 2 2 

which yields 

Re A \w\ 2 + \Tlw\ 2 + \w'\ 2 < C e 2 \w\ 2 , 

i i L 2 i i L 2 i i L 2 i i L 2 

for e and r] sufficiently small. Since \w\ 2 2 < C(|w|^ 2 + |IIw|^ 2 ) we get the estimate 
(fT2|) for e small. 

Estimate (!20l) implies the bound on the real part of the eigenvalue A. Taking the 
imaginary part of (|23l) . 

Im (X(A w,w)) = -Im (A w,d 2 FW'w) - Im (A w,dFw'} + Im (A w, dQ w). 

Hence, for rj > to be chosen, 

llmAI |w| 2 < C (\W'\ TOO |w| 2 + ?7 Iwl 2 +ri~ 1 \w'\ 2 o + r? _1 |riw;| 2 

I l l l £ 2 ' 1 'L 2 ' 1 1 L 2 ' 1 1 

Thanks to ([10]) and (JH]), we get 



|Im A| |w| 2 2 < C (e 2 + rj + ^T^ 2 ) 



2 

L2' 



Thus, choosing rj = e, we obtain the result for e small enough. □ 

3. The reduced system for the conserved integrated variables 

As stressed in the Introduction, the next step consists in estimating the conserved 
densities u in term of an appropriate weighted L 2 — norm of the conserved quantities 
z, defined by 

/X 
u(y) dy. 
-oo 

The first step is to the deduce a balance law satisfied by the variable z with source 
terms depending on Uw and w'. 
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Setting 

V -421 ^22 / V ^ ? 2 

equation (jSJ) can be rewritten as 

Att+ (AiiU + Ai2u)' = 0, 

A v + (A 21 u + A 22 v)' - q x u - q 2 v = 0. 

Let v := q 2 x q% u + v. In particular, Ci\v\ 2 < \Hw\ 2 < c 2 |-0| l2 for some C\, c 2 > 0. 
Hence, in the following, we consider the variable v in place of Uw. 
Then the couple (u, v) satisfies 

Xu + (au + A±2 v) = 0, 

(A/ r — q 2 )v+ (cu + A 22 v) + q 2 1 q 1 (au + A 12 v) = 
where 

a := A u - A X2 q 2 ~ 1 q l , c := A 21 - A 22 q 2 1 q 1 . 
With z defined in ( 1281) . we can write the above system as 
\z + az' + A\ 2 v = 0, 

{XI r -q 2 )v + (c z + A 22 v) ' + gjV ( a z' + A 12 v)' = 
Next the idea is to obtain an expression for v from the second equation and inserting 
it in the first one, in order to obtained a reduced system of viscous conservation laws 
with source terms. Since we want to derive energy estimates, it is useful to change 
variables in the first equation in order to symmetrize the term containing the first 
order derivative. 

Let do be a symmetric and positive definite matrix such that a^a is symmetric, 
as in assumption Bl. Let z := aj z. The new variable z and the variable v satisfy 

\ ~ i 1/2 / -1/2 ~v , I/ 2 A - n 

A z + a ' a [a z) + a Q ' A X2 v = (J, 

(A/ r - q 2 ) v + (c (oq 1/2 z)' + A 22 5)' + q 2 X q\ (a {a 1/2 z)' + A 12 v)' — 
Hence 

A 5 + a 5' + aj /2 A12 u = 0(1) W 5', 
(AI r -q 2 )v+ (ca" 1/2 z' + 0(1)^ + 0(l)v)' 

+q 2 1 q 1 (aa 1/2 5' + 0(l)W'z + O(l)v)' = 

where the matrix a := a\j 2 aa^ 1 is symmetric. From the second equation, using 
(fTOI). we get 

(29) {) = -(A/.-g,)" 1 ^ + ^2-^10)00 1/2 5" + 0(1)1^(5' + £5 + {)) +0(l)t)', 
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or, equivalently, using the 0(1) notation, 

(30) v = 0(l)W'(z' + ez + v) +0(l)w'. 
Plugging (I29I) in the equation satisfied by z, we get 

(31) Xz + az' -bz" = 0(l)W'(z' + ez + v) + 0(l)v' 
where 

b := Ai 2 (XI r — #2) 1 (c + q^qia) and b := a]^ 2 ba Q 1 ^ 2 . 

Estimate on z will be obtained by mutliplying (|3T|) by z and integrating. With the 
present form, we would obtain a "bad" term (z, 0(l)v'). For this reason, it is useful 
to use the relation ( 1301 to obtain the following new version of ( I3~TT) 

Xz + az' -bz" = 0(l)W'(z' + ez + v) + 0(1) (0(l)W'(z' + ez + v) +0(l)w')'. 

For A = 0, the term b represents the viscosity term given by the Chapman-Enskog 
expansion for the relaxation system. Hence, it is significant to decompose b as follows 

b = b + A61 
where matrices bo and b\ are given by 

6 := -A l2 q2 l {c + q^qia), h := A 12 (X - g 2 ) _1 ^ 1 ( c + Q^Qia)- 
Hence we get the following equation satisfied by the variable z 

Xz + az' -bo z" = Xb l z" + &. 

where b~i := a]j 2 hi a 1//2 for i = 0,1 and is appropriately defined. By assumptions 
Bl. on the reduced system, the matrix b is symmetric and positive semidefinite. 
From now on, we drop the tildas for shortness and consider the following equation 

(32) Xz + az' -b z" = Xb lZ " + e 1 + &' 2 . 

where a is a symmetric matrix, b is a symmetric, positive semidefinite matrix, and 

f 0i := 0(l)W'(z' + ez + v + w') 

( 33 ) \ s 

[ 2 := 0(l)W'(z' + ez + v) + 0{l)w' 

Lemma 3.1. Assume hypothesis C. Let X G C satisfying l[T3\) . let z be a solution 
([2J|) with a symmetric, bo symmetric, positive semidefinite and 0i, ©2 given in (E^j 
and let k be a smooth function from R™ to the set of real skew- symmetric matrices. 
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Then, for e, rj > sufficiently small, the following estimates hold: 

Re A \z\\ 2 + (z\ b z') < C{\z\l + (e + V )W\ 2 l2 + e 2 \v\\ 2 + 
Re(z',kaz'} < C(r/|Re A| \z\ 2 + s\z\ 2 w , + (r]' 1 ^ 2 + v)\ z '\\ 2 + + V' 1 \ w '\%) 

Proof. Taking the real part of the scalar product of z against (I32I) . we get 
ReA|^ 2 - -(z,daW'z) -Re(z,b z") = Re (z,Xbi z") + Re(z,G). 

having used the symmetry of a. Since 

Re (z, b z") = -Re (z', b z') - Re (z, db W'z') 

<-Re(z'Mz') + C(\z\l i + \W'\ LOO \zf L2 ). 

we obtain, thanks to ffTUj) . 

Re\\z\ 2 L2 +Re(z',b z') < C(\z\ 2 wi + e 2 \z'\\ 2 +Re(z,Xb 1 z") + Re (z, 6)) . 

The term containing b\ can be easily estimated by 

|Re (z, \b, z") | < | A| (| (z', h z')\ + \ (z, db W'z')\) 

<C\X\(\zf L2 + \z\l + \W'\ LOa \zf L2 ). 

Taking in account ( flOl) and ( fl3|) . we get 

\Re{z,Xb lZ ")\<Ce{\z'\% + \z\ 2 w ) ■ 

Therefore, we obtain 

(34) ReX\z\ 2 L2 + Re(z',b z') < C(\z\ 2 w , + e\z'\\ % + Re(z,&}). 

It remains to deal with the term with B. For what concerns 0i, using ( Tl9l) (with 
rj = 1) and ffTOl) . we have 

llte^e^l < C\z\ 2 w ,+Ce 2 {\z'\ 2 L2 + \v\ 2 l2 + \wfj 

The term with 2 can be dealt with integrating by parts 

\Re(z,Q' 2 )\ = |Re<z',9 2 )| < Ce\z\ 2 wi + C{e 2 + ^z'^ + Cs 2 \v\ 2 l2 + C^W^ 

where 7] is any positive constant and C is independent on rj. Inserting the last three 
estimates in (|34"|) . we get i. in Lemma [3. II 

Applying k to (1321) . k as defined in B2, and taking the L 2 scalar product against 
z', we get 

(35) Re (z',ka z') = -Re A (z , kz) + Re (z',kb z") + Re (z', kG) . 
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Using the eigenvalue estimate on Re A, stated in Proposition 11.51 we obtain 

Re(z',kaz') < C (^Re A||z| 2 2 + (rj'h 2 + rj) {z'^+t]' 1 |^"| 2 2 ) + Re(z', fee) 
Finally, using once more (llOl) . for the term with B there hold 

|Re (z', feGx) | < Cs\z\ 2 , + Ce 2 (\zf 2 + \v\ 2 2 + \w'\ 2 2 ) 

l ' 1 1 'w' y ' L 2 ' ' L 2 ' ' L 2 J 

|Re(z',fee' 2 )| < |Re(z",fee 2 )| + \Re (z' } 0(l)W'Q 2 }\ 

<Ce\z\ 2 +Ce 2 (\z'\ 2 2 + \v\ 2 2 )+C\w'\ 2 2 . 

W v ' ' L 2 ' ' L 2 ' 1 1 L 2 

Collecting all of these estimates, we complete the proof of Lemma [3.11 □ 

Proof of Proposition \l.b\ Choosing rj = e and summing up the estimates in Lemma 
13. 1[ we obtain, for e small enough, 

ReAN 2 2 + \z'\ 2 2 < C(\z\ 2 +e 2 \v\ 2 2 +e~ l \wf 2 ) 

L 2 ' ' L 2 v ' 'W ' ' L 2 ' ' L 2 ' 

Using (in place of a first-order Friedrichs estimate) the bound 

\v\\ 2 + W\ 2 2 <c(\Yiv\\ + W\\)<Ce 2 \u\ 2 2 

L 2 1 1 L 2 V ' L 2 1 L 2 I 1 ' L 2 



obtained in Proposition 11.51 we get 

ReAH 2 2 + M 2 < C\z\ 2 , +Ce\u\ 2 . 

L 2 1 ' L 2 ' 'W' ' ' L 2 

Hence estimate ( |T^1) holds for e small. □ 
The reduced Kawashima estimate (TH)) shows that it is possible to bound the L 2 



estimate of u in term of \z\ ,. If we are able to prove a Poincare-like inequality 
and bound the weighted norm \z\ , by small multiples of the L 2 norm of u and 



higher derivatives, we are done. This we can accomplish by changing variables in an 
appropriate way and applying a weighted energy method in the spirit of Goodman 

USE]. 

Lemma 3.2 ( [5] ). Let a = a(W) and b = b(W) be symmetric matrices, b(W) > ; 
with one eigenvalue cto of a close to zero and the others strictly negative or positive 
(and uniformly separated from a^). Then, there exist smooth, real matrix-valued 
functions r = r(x), I = £(x), £(x) r(x) = I for any x, satisfying, for some C,c> 0, 

(£r') pp = 0, \£'\,\r'\ <C\W'\; 



(36) Ear = diag (a_, a p , a+) 
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with a p scalar, a_, a + symmetric square matrices (with dimensions p — 1 and n — p 
respectively), at- < — c < < c < a + ; and 

(37) Re£br>-Cs. 

Proof. Since a is symmetric, it is possible to find an othonormal transformation 
uj = uo{W) such that u 1 auo is (block-) diagonal with the decomposition given in the 
righthand side of ([55]) . The spectral separation assumption guarantees the posi- 
tivity /negativity of ct + /a_. Moreover, the matrix ^bu is positive semidefinite, 
Re6 > 0. 

Let Up denote the pth column of to and 7 = 7(0;) be the solution of the first order 
linear differential equation 

(38) i = -(u p -du Jp W') 1 , 7 (0) = 1, 
or, equivalently, set 

(39) 7(2) := exp (J u* p (W) dcu p (W) W dy 
Define the matrix r and i as 

(40) r{x) := u{W) diag 7 (x), J n _ p ), t{x) := r" 1 ^). 

Clearly estimates on |r'| and hold and 

£ar = diag (J p _i, 7" 1 , J n _ p )diag (a_, a p , a+)diag (J p _i, 7, J n __ p ) 

= diag (a_, a p , a+), 
hence £ and r still block-diagonalize a in the manner claimed. Moreover 

(£r')pp = j' 1 ^ ■ ("fujp)' = u p ■ (du p W' - (ujp ■ dup W')uj p ) = 0, 

since ui p has norm equal to 1. 

By (I3"9l . it follows that 7 = 1 + 0(e); hence bound (jBTjl follows from assumptions 
on 6 and continuity. □ 

Proof of Proposition \ 1 . 7[ To prove (jl5p . it is sufficient to establish the corresponding 
result for ( := £ z with £ given in Lemma l3~2l Left multiplying (|32|) by we get the 
eigenvalue equation for ( 

(41) AC + o^r'C + = 
where 

a:=£ar, (3 := £b r, E := Hi + E' 2 
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and 

~ i: =0(l)\W'\{e( + C + v' + w'), 

E 2 :=0(l)\W'\(t + v) +0(l)eC + 0(l)w'. 
Following [3j, set po(x) : = 1 for any x, and define the two weights p± as the solutions 
to the Cauchy problem 

(42) p ' ± = T M\W'\c- 1 p ± , p±(0):=l, 

where c is given in Lemma 13.21 and M is a constant to be chosen later. Therefore, 
for e so small that O(Me) < 1, 

(43) p±(x) =exp(±^ M\W'(£)\c- x dZ) = 1 + O \M J^\W' d£) =0(1), 
and 

(44) p' j (x) = 0(l)\W'(x)\, 3 e{-,0,+}. 
Let p = p(x) be the block diagonal matrix defined by 

p(x) := diag(p_(x) I h , p (x), p+(x) h) 

where I n denotes the identity n x n matrix. Taking the real part of the L 2 — scalar 
product of p( against ( 14T|) . we get 

(45) Re A (p C, C) + Re(p(,a£ r'() + Re (p (, a (') - Re (p (, (3 (") = Re(p(,~). 

The weights po, p± are positive and 0(1), hence (pCX) 1 ^ 2 is equivalent to \(\ l2 - 
Both p and pa are symmetric, hence 

Re (p C, a O = Re (C, p a (') = ~Re (C, (p' a + p da W) (). 

By (J9]) and (TTTj) . we have the key fact^ 

da W < -C\W'\ 

for some C > 0. By definition of p±, we have also 

p'± a± + P± da± W = TP± (M I W"|c- V - da± W). 

Thus, for M sufficiently large, there exists C > 0, independent on e, such that 
p' a + pdaW' = diag (p'_ ct_ + p_ da_ W 7 ', da W, p' + a + + p + da + W') 
< -C\W'\diag(M,l } M). 

3 Indeed, this is what drives the Goodman estimate; see [31 3]. 
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Decomposing ( as (C- , Co? C+) an d setting ( := (C-?C+)> we get the "good" term 

Re(p(,aC)>C [ (M\(\ 2 + \( \ 2 )\W'\dx. 
Jr 

Next, let us deal with the "bad" term (p(,a£r'()- Since (£r') pp = 0, there holds 

\Re(p(,a£r'()\<C [ \(\ 2 \W'\dx 



Hence, by choosing M large enough, we get from (1451) 

ReA IC| 2 2 + IC -Re(K,/?C"> < C|Re(pC,S}| 

Since 

Re (p C, (3 C"> = -Re (p P O ~ Re <P C, C) - Re <p' C, /? 
the term with /3 can be estimated by 

Re(pC,/50<C £ (|C| 2 w , + |C'| 2 L2 ) 
having used ( J37l) . Hence, we obtain 

(46) ReA|C| 2 L2 + lC <^(Kt + ICt 2 )+C|Re(pC,S)| 

Given r\ > 0, recalling ( fTOl) . we deduce 

|He(pC,Hi>| < C^ + ^KI^ +Cr / -V(|C'| 2 L2 + \v\% + [wfj 

with C independent on rj. For what concerns the term with S 2 , integrating by parts 
and using ( pHl) . there holds 

|Re(pC,H' 2 >| = |Re(p'C,H 2 )| + |Re(pC',2 2 )| < |Re(0(l)|W| C,H 2 )| + |Re(pC',H 2 )| 

For any rj > 0, estimating one by one the terms in 5 2 , we obtain 

|Re (0(1)|W'| C, H 2 }| < C(e + ??)|C| 2 „ + Ce 2 (|C'| 2 2 + \v\ 2 J + CVVK| 2 2 , 

| Re (pC, S 2 }| < C77|Ci 2 w; + C(e + tT^ 2 + ^ICI^ + C^e 2 \v\\ % + CV V| 2 2 - 
Choosing rj = e and summing up, we get 

|Re (pC,H 2 )| < Ce(\C\l + |C'| 2 2 + \v\%) + Ce^w'^. 
Inserting these estimates in fH6|) . we get, for £ sufficiently small, 

ReA|C|^ 2 + \C\ 2 WI < C(\C\ 2 L2 +e\v\% +^V| 2 l2 ). 
Since ( = £ z, from the above estimate we deduce 

ReA|< 2 + \z\ 2 <C(e\z'\ 2 + e\v\ 2 + e~W ). 
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for e sufficiently small. Recalling that z' = u, estimate ( 115}) is proved. 



□ 



Appendix A. Structural hypotheses 

In this Appendix, we briefly discuss the structural hypotheses of the introduction, 
verifying the assertions of Remark 1 1.21 that A1-A2 and B1-B2 follow from conditions 
(i)-(ii) of the remark (i.e., partial simultaneous symmetrizability plus genuine cou- 
pling) together with the assumed structure Q = (0 n ,q). 

Lemma A.l. Let Q = (0 n ,q). Then, (i)-(ii) of RmkEIM imply A1-A2 and B1-B2. 

Proof. These follow by more general results of Yong We give a proof for com- 

pleteness. As all properties are coordinate-independent properties of the lineariza- 
tion about constant states, we may without loss of generality take A block-diagonal. 
For, TA° is block-lower triangular for T block-upper triangular, whence TA°T* is 
symmetric block- diagonal, and a left symmetrizer for the system obtained by the 
change of coordinates w -> (T*)~V A -> (T*)' l AT*, Q -> (T*)- l QT*. 

Observing that A := (A )" 1 is a right symmetrizer if A is a left symmetrizer, 
we obtain 

A°w t + Aw x = Qw, 

where A is symmetric positive definite and block-diagonal, A is symmetric, and 
q ~^ symmetric with q < 0. (Note: the latter key fact follows by 

dQ = (T*)~M<5 T*A°, 

the fact that T* is block-lower triangular, and that the first block row of dQ by 
assumption vanishes.) Rewriting, we have 

w t + Aw x = Qw, 

where 

\(A° 22 )-'A 21 (A^AJ' ^ vO (A^qJ- 
In these coordinates, one readily computes that 

a = (A^-'An, b = -(A° n )^A* 12 qA* 12 , 

hence ker6 = kerA^, and genuine coupling for the reduced system is the condition 
that no eigenvector of oq = (A® 1 )~ 1 Au lie in kerAi2, the same condition as for 



Symmetrizability is not explicitly stated in [10] . but is clear from the development. 
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genuine coupling of the full system, and is a left symmetrizer for the reduced 
system with A^bo = — Al 2 qAl 2 symmetric positive semidefinite since q is symmetric 
negative definite. □ 
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